Charge Ordering in Half-Doped Manganites: Weak Charge Dis- 
proportion and Leading Mechanisms 



Dmitri Volja 1 - 2 , Wei-Guo Yin 1 ^ and Wei Ku 1 ' 2 ^ 



1 Condensed Matter Physics and Materials Science Department, Brookhaven National Laboratory, Upton, New York 
11973, USA 

c"^ ^ Physics Department, State University of New York, Stony Brook, New York 11790, USA 



o 



o- 
o 



PACS 75 . 47 . Lx - Manganites (magnetotransport materials) 

PACS 71.45.Lr - Charge-density waves - collective excitations 

PACS 71. 10. Fd - Lattice fermion models (Hubbard model, etc.) 

PACS 71 .30.+h - Metal-insulator transitions and other electronic transitions 

Abstract. - The apparent contradiction between the recently observed weak charge disproportion 
and the traditional Mn 3+ /Mn 4+ picture of the charge-orbital orders in half-doped manganites is 
resolved by a novel Wannier states analysis of the LDA+C/ electronic structure. Strong electron 
itinerancy in this charge-transfer system significantly delocalizes the occupied low-energy "Mn 3+ " 
Wannier states such that charge leaks into the "Mn 4+ "-sites. Furthermore, the leading mechanisms 
of the charge order are quantified via our first-principles derivation of the low-energy effective 
Hamiltonian. The electron-electron interaction is found to play a role as important as the electron- 
lattice interaction. 
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Introduction. — The exploration of interplay among 
distinct orders lies in the heart of condensed matter 
physics and materials science, as this interplay often give 
rises to tunable properties of practical applications, such 
as exotic states and colossal responses to external stimuli. 
Manganese oxides such as Lai-^Ca^MnOs, which host 
rich charge, orbital, spin, and lattice degrees of freedom, 
have thus attracted great attention [T|. In particular, the 
vastly interesting colossal magnetoresistance (CMR) ef- 
fect for x ~ 0.2 — 0.4 exemplifies rich physics originating 
from proximity of competing orders. In a slightly more 
doped system (x — 0.5), all these orders coexist in an in- 
sulating state 2,, providing a unique opportunity for a 
clean investigation of the strength and origin of each or- 
der. Therefore, the study of half-doped manganites is key 
to a realistic understanding of the physics of manganites 
in general and the CMR effect in particular [3j|4]. 

The peculiar multiple orders in half-doped manganites 
have long been understood in the Goodenough model [2] 
of a Mn 3+/4+ checkerboard charge order (CO) with the 
occupied Mn 3+ e g orbitals zigzag ordered in the CE-type 
antiferromagnetic background [2J. Pertaining to CMR, it 
is broadly believed that a key is the emerging of nanoscale 
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charge-ordered insulating regions of Goodenough type at 
intermediate temperatures, which could melt rapidly in 
the magnetic field [H[S]. 

Nonetheless, the simple yet profound Goodenough pic- 
ture has been vigorously challenged thanks to recent ex- 
perimental observations of nearly indiscernible charge dis- 
proportion (CD) in a number of half-doped manganites 
[5Hl3j. Such weak CD has also been observed in first- 
principles computations [141115) and charge transfer be- 
tween Mn and O sites was reported as well [TS]. In 
essence, these findings have revived a broader discussion 
on the substantial mismatch of valence and charge in most 
charge-transfer insulators. Indeed, extensive experimen- 
tal and theoretical effort has been made in light of the 
novel Zener-polaron model [M53"] in which all the Mn sites 
become equivalent with valence being +3.5. Amazingly, 
most of these investigations concluded in favor of two 
distinct Mn sites as predicted in the Goodenough model 
jTTHI31[T51|T7] , calling for understanding the emergence of 
weak CD within the 3+/4+ valence picture. 

Another closely related crucial issue is the roles of dif- 
ferent microscopic interactions in the observed charge- 
orbital orders, in particular the relevance of electron- 
electron (e-e) interactions in comparison with the well- 
accepted electron-lattice (e-l) interactions. For example, 
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An, the difference in the electron occupation number b( 
tween Mn 3+ and Mn 4+ states, was shown to be small in a 
e-e only picture |18) ; this is however insufficient to explai 
the observed weak CD, since the established e-l interae 
tions will cause a large An [T5]. Moreover, despite tb 
common belief that e-l interactions dominate the gener; 
physics of e g electrons in the manganites [HHU], a recer 
theoretical study showed that e-e interaction plays a 
essential and leading role in ordering the e g orbitals in tb 
parent compound. It is thus important to quantify tb 
leading mechanisms in the doped case and uncover the e: 
fects of the additional charge degree of freedom in genera 
In this Letter, we present a general, simple, yet quant 
tative picture of doped holes in strongly correlated charg< 
transfer systems, and apply it to resolve the above cor 
temporary fundamental issues concerning the charge orde 
in half-doped manganites. Based on recently develope 
first-principles Wannier states (WSs) analysis 21,22,34 
the LDA+J7 electronic structure in prototypical Ca-dop< 
manganites, the doped holes are found to reside primari 
in the oxygen atoms. They are entirely coherent in sho 
range [23], forming a Wannier orbital of Mn e g symm 
try at low-energy scale. This hybrid orbital, together wil 
the unoccupied Mn 3c? orbital, forms the effective "Mn e, 
basis in the low-energy theory with conventional 3+/4 
valence picture, but simultaneously results in a weak C 
owing to the similar degree of mixing with the intrins 
Mn orbitals, thus reconciling the current conceptual co: 
tradictions. Moreover, our first-principles derivation of tl 
low-energy interacting Hamiltonian reveals a surprising 
essential role of e-e interactions in the observed charge o 
der, contrary to the current lore. Our theoretical methc 
and the resulting simple picture provide a general fram 
work to utilize the powerful valence picture even with we; 
CD, and can be directly applied to a wide range of dop( 
charge-transfer insulators. 

Small CD vs. 3+/4+ valence Picture. To pr 

ceed with our WSs analysis, the first-principles electron 
structure needs to reproduce all the relevant experiment 
observations, including a band gap of ~ 1.3 eV, CE-ty] 
magnetic and orbital orders, and weak CD, as well as tv 
distinct Mn sites. We find that the criteria are met by tl 
LDA+U (8 eV) [HIEM] band structure of the prototypic 
half- hoped manganite La 1 / 2 Ca 1 / 2 Mn03 based on the r 
alistic crystal structure [35] supplemented with assume 
alternating La/Ca order. Hence, a proper analysis of tb 
LDA+U electronic structure is expected to illustrate tl 
unified picture of weak CD with the Mn 3+ /Mn 4+ assig: 
ment, which can be easily extended to other cases. In 
practice, we shall focus on the most relevant low-energy 
(near the Fermi level Ep) bands — they are 16 "e g " spin- 
majority bands (corresponding to 8 "spin-up" Mn atoms 
in our unit cell) spanning an energy window of 3.2 eV, 
as clearly shown in Fig.QJa). For short notation, the Mn 
bridge- and corner-sites in the zigzag ferromagnetic chain 
are abbreviated to B- and C-sites, respectively. 




Fig. 1: (Color online) (a) LDA+[/ Band structures (dots). The 
(red) lines result from the Wannier states analysis of the four 
occupied spin-majority Mn 3d-derived bands, (b) An occupied 
B-site ("Mn 3+ ") Wannier orbital in a spin- up (up arrow) zig- 
zag chain, showing remarkable derealization to the neighbor- 
ing Mn C-sites. (c) Low-energy Mn atomic-like Wannier states 
containing in their tails the integrated out O 2p orbitals. 



The simplest yet realistic picture of the CO can be ob- 
tained by constructing occupation-resolved WSs (ORWSs) 
from the four fully occupied bands, each centered at one 
B-site as illustrated in Fig. QTa). This occupied B-site e g 
Wannier orbital of 3x 2 — r 2 or 3y 2 — r 2 symmetry (so for- 
mal valency is 3+) contains in its tail the integrated out 
O 2p orbitals with considerable weight, indicative of the 



p-2 



Charge Ordering in Half-Doped Manganites 



charge-transfer nature [IS]- Moreover, this "molecular or- 
bital in the crystal" extends significantly to neighboring 
C-sites on the same zig-zag chain. Therefore, although by 
construction the C-site e g ORWSs (not shown) are com- 
pletely unoccupied (so formal valency is 4+), appreciable 
charge is still accumulated within the C-site Mn atomic 
sphere owing to the large tails of the two occupied OR- 
WSs centered at the two neighboring B-sites. Integrating 
the charges within the atomic spheres around the B- and 
C-site Mn atoms leads to a CD of mere 0.14 e, in agree- 
ment with experimental 0.1 — 0.2 e [7HT2"]. In this simple 
picture, one finds a large difference in the occupation num- 
bers of the ORWSs at the B- and C- sites (An — I), but 
a small difference in real charge. That is, the convenient 
3+/4+ picture is perfectly applicable and it allows weak 
CD, as long as the itinerant nature of manganites is incor- 
porated via low-energy WSs rather than standard "atomic 
states." 

In comparison, to make connection with the conven- 
tional atomic picture and to formulate the spontaneous 
symmetry breaking with a symmetric starting point (c.f. 
the next section), we construct from the 16 low-energy 
bands "atomic-like" WSs (AWSs) of Mn d 3z 2„ r 2 and 
d x 2_ y 2 symmetry, as shown in Fig. [ljb). In this pic- 
ture, both B- and C-site AWSs are partially occupied with 
An = 0.6. Now weak CD results from large hybridization 
with O 2p orbitals, which significantly decreases the charge 
within the Mn atomic sphere. Obviously, this picture is 
less convenient for an intuitive and quantitative under- 
standing of weak CD compatible with the 3+/4+ picture 
than the previous one, as the latter builds the information 
of the Hamiltonian and the resulting reduced density ma- 
trix into the basis. On the other hand, it indeed implies 
that strong charge-transfer in the system renders it highly 
inappropriate to associate CD with the difference in the 
occupation numbers of atomic-like states. 

Furthermore, we find that the above conclusions are 
generic in manganites by also looking into the two end lim- 
its of Lai^Ca^MnOs (x = 0, 1). As shown in Fig.0 Mn d 
charge (within the Mn atomic sphere) is found to change 
only insignificantly upon doping, in agreement with ex- 
periments [?1[26]. This indicates that doped holes reside 
primarily in the oxygen atoms, but are entirely coherent 
in short-range and form additional effective "Mn e ff " or- 
bitals in order to gain the most kinetic energy from the 
d-p hybridization, as shown in Fig. [IJa)-(b), in spirit sim- 
ilar to hole-doped cuprates [23 . This justifies the present 
simplest description of the charge-orbital orders with only 
the above Mn-centering e g WSs [2"T] . 

Leading Mechanisms. — To identify the leading 
mechanisms of the charge- orbital orders in a rigorous for- 
malism, we proceed to derive a realistic effective low- 
energy Hamiltonian, H cS , following our recently devel- 
oped first-principles WS approach [3TJ|33]. As clearly 
shown in Fig. [TJ the low-energy physics concerning charge 
and orbital orders is mainly the physics of one zig-zag FM 
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Fig. 2: The calculated d charge within the Mn atomic sphere 
for Lai-zCa^MnOs (x = 0, 1/2, and 1). The results are shown 
in terms of (i) B-site (filled symbols), C-site (empty symbols), 
and site-average (half-filled symbols); (ii) spin- majority (up tri- 
angles), spin-minority (down triangles), and total charge (dia- 
monds). The lines are guides to eye. 

chain, since electron hopping between the antiferromag- 
netically arranged chains is strongly suppressed by the 
double-exchange effect [IH1I2E1I22- Our unbiased first- 
principles analysis of the 16-band one-particle LDA+U 
Hamiltonian in the above AWS representation reveals 

HSU2SH22U3BI 

<ij)77' i 

+ U eS ^ n n n ii + v zl ni "J 

i <ij> 

- SZX^ niQu + TfQx + T.zQ^) (1) 

in addition to the elastic energy K{{Qi}). Here d| 7 and 
dij are electron creation and annihilation operators at site 
i with "pseudo-spin" 7 defined as | f) = |3z 2 — r 2 ) and 
I D — \y 2 ~ x 2 ) AWSs, corresponding to pseudo-spin oper- 
ator T? = (4r d u + d u d it)/ 2 and T i = (4t d it- d h d n.)/ 2 - 
Ui = °4t^it dlidu is the electron occupation number. 
The in-plane hoppings are basically symmetry related: 
t$ = t/4, tfj- = 3t/4, t^ = t$ = ±V3t/4 where the signs 
depend on hopping along the x or y direction. E z stands 
for the oxygen octahedral-tilting induced crystal field. U e g 
and V are effective on-site and nearest-neighbor e-e in- 
teractions, respectively, g is the e-l coupling constant. 
Qi = (Qii: Q21, Q31) is the standard octahedral-distortion 
vector, where Qu is the breathing mode (BM), and Q21 
and Q 3i are the Jahn- Teller (JT) modes [Mll9l[28l[29ll36] . 
In Eq. (1) the electron- lattice couplings have been con- 
strained to be invariant under the transformation of the 
cubic group |36j . 

The effective Hamiltonian are determined by match- 
ing its self-consistent Hartree-Fock (HF) expression with 
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Table 1: Contributions of different terms to the energy gain 
(in units of meV per Mn) due to the CO formation in self- 
consistent mean-field theory. BM (JT) denotes the contribu- 
tion from electronic coupling to the breathing (Jahn- Teller) 
mode. K denotes that from the elastic energy. 

Qi Total U cS V t BM JT K 
-13 -11 -15 12 

realistic -127 -22 -42 71 -42 -113 24 



jjLDA+u [2T1[34] ow i ng to the analytical structure of the 
LDA+C/ approximation [30] . An excellent mapping re- 
sults from t = 0.6 eV, E z = -0.08 eV, U cS = 1.65 eV, 
V = 0.44 eV, and g = 2.35 eV/A. These numbers are 
close to those obtained for undoped LaMnOs [2T] (exclud- 
ing V, which is inert in the undoped case), and indicates 
that the spin-majority e g electrons in the manganites are 
still in the intermediate e-e interaction regime with compa- 
rable e-l interaction. Note that U c s should be understood 
as an effective repulsion of corresponding Mn-centered WS 
playing the role of "d" states, rather than the "bare" d-d 
interaction U — 8 eV [21] . Furthermore, the rigidity of the 
low-energy parameters upon significant (a; = 0.5) doping 
verifies the validity of using a single set of parameters for 
a wide range of doping levels, a common practice that is 
not a priori justified for low-energy effective Hamiltoni- 
ans. Clearly, the observed optical gap energy scale of ~ 2 
eV originates mainly from U c s instead of the JT splitting 
widely assumed in existing theories [51 [211] . 

Now based on the AWSs, CO is measured by An = 
(nigs) — ("jec) = 0.6. It deviates from unity because the 
kinetic energy ensures that the ground state is a hybrid of 
both B- and C-site AWSs, like in the usual tight-binding 
modeling based on atomic orbitals. However, since the 
AWSs considerably extend to neighboring oxygen atoms, 
the actual CD is much smaller than An. 

With the successful derivation of H eS , the microscopic 
mechanisms of the charge-orbital orders emerge. First of 
all, note that the kinetic term alone is able to produce the 
orbital ordered insulating phase [T81I29] : Since the intersite 
intcrorbital hoppings of the Mn e g electrons along the x 
and y directions have opposite signs, the occupied bonding 
state is gapped from the unoccupied nonbonding and an- 
tibonding states (by t and 2t, respectively) in the enlarged 
unit cell. As for orbital ordering, the Mn d y 2_ z 2 (d x 2_ z 2) 
orbital on any B-site bridging two C-sites along the x (y) 
direction is irrelevant, as the hopping integrals involving 
it is vanishing. That is, only the d 3x 2_ r 2 (d 3y 2_ r 2) orbital 
on that B-site is active and the B-sites on the zigzag FM 
chains have to form an "ordered" pattern of alternating 
(3x 2 — r 2 )/(3y 2 — r 2 ) orbitals. However, the kinetic term 
alone give only An = 0. Clearly, CO is induced by the 
interactions, U c s, V, or g. 

To quantify their relative importance for CO, we calcu- 
late their individual contributions to the total energy gain 



with respect to the aforementioned An = but orbital- 
ordered insulating state in the self-consistent mean-field 
theory. The results are listed in Table [TJ The first row 
obtained without lattice distortions provides a measure of 
the purely electronic mechanisms for CO. Interestingly, al- 
though the tendency is weak (—13 meV), e-e interactions 
all together are sufficient to induce a CO, consistent with 
the results of our first-principles calculations. The second 
row is obtained for the realistic lattice distortions, which 
shows a dramatic enhancement of CO by the JT distor- 
tions (—113 meV), given that only half of the Mn atoms 
are JT active. Together with a —42 meV gain from the 
BM distortion, the e-l interactions overwhelm the cost of 
the kinetic (71 meV) and elastic (24 meV) energy by —63 
meV, further stabilizing the observed CO. Nevertheless, 
the —64 meV energy gain from the overall e-e interac- 
tions accounts for half of the total energy gain, illustrating 
clearly their importance to the realization of the resulting 
An = 0.6. Indeed, a further analysis reveals that e-l cou- 
plings alone (U e g = V = 0) produce An ~ 0.3, only half 
of the realistic An, manifesting the necessity of including 
e-e interactions. 

Further insights can be obtained by considering the in- 
dividual microscopic roles of these interactions acting to 
the kinetic-only starting point. First, infinitesimal U e g or 
g can induce CO, as a result of exploiting the fact that 
the B- (C-) site has one (two) active e g orbital: (i) [7 q 
has no effect on the B-sites; therefore, U e g pushes the e g 
electrons to the B-sites, in order to lower the Coulomb 
energy on the C-sites [18:. This is opposite to its normal 
behavior of favoring charge homogeneity in systems such 
as straight FM chains (realized in the C-type antiferro- 
magnet). (ii) It is favorable to cooperatively induce the 
(3x 2 — r 2 )/(3?/ 2 — r 2 )-type JT distortions on the B-sites 
and the BM distortions on the C-sites in order to mini- 
mize elastic energy [2j . These lattice distortions lower the 
relative potential energy of the active orbitals in the B- 
sites, also driving the e g electrons there. Hence, U e Q and 
e-l interactions work cooperatively in the CO formation. 

Unexpectedly, we find that V alone must be larger than 
V c ~ 0.8 eV to induce CO. This is surprising in compar- 
ison to the well-known V c = for straight FM chains. 
The existence of V c is in fact a general phenomenon in 
a "pre-gapped" system. Generally speaking, in a sys- 
tem with a charge gap, A , before CO takes place (e.g. 
the zigzag chain discussed here), forming CO always costs 
non-negligible kinetic energy due to the mixing of states 
across the gap. As a consequence, unlike the first order 
energy gain from U e g and g, the second order energy gain 
from V is insufficient to overcome this cost until V is large 
enough (of order A ). In this specific case, V = 0.44 eV is 
insufficient to induce CO by itself, but it does contribute 
significantly to the total energy gain once CO is triggered 
by U c s or g, as discussed above. 

It is worth mentioning that the contribution from the 
E Z T Z term is neglected from Table 1 because of the small 
coefficient of E z ~ 8 meV, consistent with the previous 
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study for undoped LaaMnOs [21]. In perovskites, the tilt- 
ing of the oxygen octahedra could yield the Jahn- Teller- 
like distortion of GdFcC>3 type, which can be mathemat- 
ically described by the E Z T Z or E X T X term. However, in 
the perovskite manganites these terms are shown here and 
in Ref. [2T] to be negligibly small. In addition, in half- 
doped manganites, the pattern of the orbital order is pre- 
dominantly pined by the zig-zag pathway of the itinerant 
electrons and thus the effect of tilting is less relevant. On 
the other hand, the E Z T Z is very effective to explain the 
zig-zag orbital ordering of (x 2 — z 2 )/(y 2 — z 2 ) pattern in 
half-doped layered manganites, such as Lao.sSri.5Mn04 
Ref. [35 1, where pseudo-spin-up (the 3z 2 — r 2 orbital) is 
favored by a much stronger E z due to the elongation of 
the oxygen octahedral along the c-axis. Described via 
the pseudo-spin angel, 9 — arctan(T x /T z ) [21] . E z signif- 
icantly reduces from ±120° [i.e., (3x 2 - r 2 )/(3y 2 - r 2 )] 
to ±30° [i.e., (x 2 - z 2 )/(y 2 - z 2 ). 

The present results would impose stringent constraints 
on the general understanding of the manganites. For ex- 
ample, Zener polarons were shown to coexist with the CE 
phase within a purely electronic modeling of near half- 
doped manganites [3T]. However, to predict the realistic 
phase diagram of the manganites, one must take into ac- 
count the lattice degree of freedom. In fact, when propos- 
ing the CE phase, Goodenough considered its advantage 
of minimizing the strain energy cost. Note that the pre- 
vious HF theory indicated a decrease of the total energy 
by 0.5 eV per unit cell with Zener-polaron-like displace- 
ment [13]. The present LDA+U calculations reveal an 
increase of the total energy by 1.07 eV per unit cell with 
the same Zener-polaron-like displacement. This discrep- 
ancy is quite understandable from the characteristics of 
the LDA functional, which favors covalent bond, while the 
HF approximation tends to over localize the orbital and 
disfavor bonding. To resolve the competition between the 
CE phase and the Zener polaron phase, real structural 
optimization is necessary and will be presented elsewhere. 
As another example, in the absence of e-l interactions the 
holes were predicted to localize in the B-site-type region 
(i.e., the straight segment portion of the zigzag FM chain) 
in the C^Ei-^ phase of doped E-type manganites (32j . 
However, since the B-sites are susceptible to the JT dis- 
tortion, they are more likely to favor electron localization 
instead; future experimental verification is desirable. 

Summary. — A general first-principles Wannier func- 
tion based method and the resulting valence picture of 
doped holes in strongly correlated charge-transfer systems 
are presented. Application to the charge order in half- 
doped manganites reconciles the current fundamental con- 
tradictions between the traditional 3+/4+ valence picture 
and the recently observed small charge disproportion. In 
essence, while the doped holes primarily resides in the oxy- 
gen atoms, the local orbital are entirely coherent following 
the symmetry of Mn e g orbital, giving rise to an effective 
valence picture with weak CD. Furthermore, our first- 



principles derivation of realistic low-energy Hamiltonian 
reveals a surprisingly important role of electron-electron 
interactions in ordering charges, contrary to current lore. 
Our theoretical method and the resulting flexible valence 
picture can be applied to a wide range of doped charge- 
transfer insulators for realistic investigations and interpre- 
tations of the rich properties of the doped holes. 
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